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Abstract 

We consider infinite particle system on the positive half-line moving independently of each other. When 
a particle hits the boundary it immediately disappears, and the boundary moves to the right on some fixed 
quantity (particle size). We study the speed of the boundary movement (growth). Possible applications - 
dynamics of the traffic jam growth, growth of thrombus, epitaxy. Nontrivial mathematics is related to the 
correlation between particle dynamics and boundary growth. 

1 Introduction 

Hitting the boundary by the random walk is a classical problem of probability theory. More complicated is the 
problem of hitting the moving boundary, see pQ. Here we consider the case when the movement of the boundary 
is correlated with the movement of the particles. Namely, we consider infinite number of particles, performing 
random walks on the positive half-line, which stick to the boundary when hit it. Moreover, when a particle hits 
the boundary, the boundary moves to the right on some quantity, "size" of the particle. That is the boundary 
moves due to accumulation of particles on it. As far as we know, such problems were not considered earlier. 

Similar problems one can often encounter in applications - growth of thrombus, epitaxy and other methods 
of coating the surface with metal particles. In many cases the effect of correlation between the movements of 
particles and the boundary can be neglected. Here, on the contrary, we consider the influence of this correlation 
on the growth speed of the boundary. Such effects can be observed while quick formation of traffic jams. 

In the paper we present two results. First one for nonzero drift of the particles, second - for zero drift. In 
the first case one can find exact asymptotic rate of growth, in the second case one can get only the order of 
growth. 

2 Formulation of the problem and the results 

We consider infinite system of particles on the half-line, which, while hitting the boundary, adhere to it (accu- 
mulate on it). As a result the boundary moves to the right. Now we give exact formulation. 
Assume that at time on the half- line i? + the random configuration of particles at the points 

< xi = xi(0) < ... < x n = x n (0) < ... 

is distributed as the point Poisson process with density A. 

Each particle, before hitting the boundary, moves independently of other particles as 

Xi(t) = Xj — Vt + Wj(t) (1) 

where v > is the constant drift, and Wi(t) is the standard Wiener process with zero mean. At the hitting 
moment the particles are absorbed by the boundary and disappear, moreover the boundary moves to the right 
on some 5, the "particle size". 

More exactly, the boundary which at time t = coincides with point G R+, also moves. The law £(i) of 
its movement models the sticking of particles to the boundary. Namely, £(0) = 0, and £(i) is a non-decreasing 
point-wise constant function with jumps 



i 



at random time moments 

< ti < ... < tj < ... 

where 6 > is some constant, the "particle size", the moments tj and the non-negative integers kj being defined 
by the following recurrent procedure. Let t\ be the first instant, when one of the particles hits 0, then k\ is the 
minimal integer k > such that in the interval (0, kS] at this moment t\ there are not more than k — 1 particles. 
For example, k = 1 if and only if at time t\ there are no particles inside the interval (0,(5]. Further on by 
induction, t„+i is defined as the first moment when one of the remaining particles hits the point (fci + ... + k n )S, 
and k n+ \ is defined correspondingly as the minimal integer k > such that inside the interval (£(t n ),k5] at 
time t n+ i there is no more than fc — 1 particles. 
Below the following statement will be proved. 

Lemma 1 If A < S^ 1 , then with probability 1 all fcj and ti are finite, and moreover ti oo. On the contrary, 
if A > 5~ x , the boundary reaches infinity for finite time with positive probability. 

Now we can formulate the main result. Let N = N(t) be the number of particles, absorbed by the boundary 
to time t, then £ = = SN(t) is the coordinate of the particle at time t. 

Theorem 2 Let A < <5 _1 . Ifv>0, the particle movement is asymptotically linear, that is as t — > oo a. s. 



If v = 0, then for t sufficiently large 
for some constants 



t 1-SX 



C <^1<C 
yt 



Ay - < C < C < oo 



Remark 3 We will show later that under the condition that the boundary does not move at all and v = 0, 
then the mean number of particles hitting the boundary during time t equals Ay^. From the second assertion 
of the theorem it follows in particular that the mean number of particles hitting the boundary is asymptotically 
greater. We conjecture that MN ^ tends to some constant as t — > oo, and moreover the distribution of the 



random variable ^fi- tends to some continuous distribution as t — > oo. 



Vt 



3 Proofs 

Proof of Lemma [J For a given configuration X = {xi(0)} denote p(X,n) the probability that the first 
particle hitting the boundary is the particle x n (0). 

To prove the first part of the Lemma introduce the following auxiliary model. At time t = for a given 
configuration on the positive axis we add Poisson configuration on the negative axis Y, having the same density 
A. Then we will get Poisson configuration on all line R, all these particles have red color. Moreover, there is 
one more (blue) particle, which we put with probability p(X, n) at the point x n (0) > 0. The particles move 
according to ([1}, independently of each other. At time t = the boundary is at the point £'(0) = 0. The 
random variable t\ is defined as the first moment when the blue particle hits 0, and k[ are defined as earlier. 
The difference from the basic model is that the boundary does not influence other particles at all. Thus at any 
moment t the red particles g/j(t), — oo < j < oo, have Poisson field distribution (Doob-Dobrushin theorem, [3]). 
Moreover at random time t[ the red particles have Poisson field distribution. 

Let £'(£) be the position of the boundary at time t. Then for all t a. s. £'(t) > £(t). In particular, k[ > k\ 
a. s. In fact, let 51 be the probability space of the auxiliary model. The basic model is defined on the same 
probability space, as it can be obtained from the auxiliary model by deleting the glue particle and the particles 
situated to the left of 0. Hence in the auxiliary model the number of collisions with the boundary is not less 
than in the basic model. 



2 



Let us prove that k[ is finite a. s. for any n. The random variable k[ is defined by the stopping time of the 
following Markov chain mi, m„ ... with discrete time and state space Z + . The state is an absorbing state. 
Denote mithe number of particles in the interval [0, 6] at time t%, and mi the number of particles in the interval 
[5, (mi + 1)5], if mi > 0, m.3 the number of particles in the interval [(mi + 1)6, (mi + 1)<5 + m,2<5], if mi > 0, 
m-2 > 0, and so on. The sequence m, constitute time homogeneous Markov chain with transition probabilities 

p mi m i+1 - TOi+i , 

The boundary stops when mi = 0. One step increment is 

M(mj+i — mi\rrii = k) = k(X5 — 1) 

The chain hits with probability 1, if X5 — 1 < 0. 

By induction, using similar trick with appending a new blue particle at any moments t[ , the finiteness of all 
k[ follows. 

The second statement of the Lemma is not used in the proof of the Theorem, and it will be convenient to 
prove it at the end of the paper after the proof of the Theorem. 



Idea of the proof of the Theorem One can write down two equalities relating random variables £(t) and 
N(t): 

^(t) = N(t)S (2) 

N(t)=n (t,m) + ni(t^(t)) (3) 

where no(t, £) is the number of particles which were at time inside the interval [0, £] and were absorbed by 
the boundary before time t, ni(i, £) is the number of particles which were at time to the right of the point £ 
and were absorbed by the boundary before time t. We can write the following equation 

It will appear useful if we shall have good expressions for no and ri\ . 

The first result of the theorem can be explained for the case when for nonzero drift there is no fluctuations. 
More exactly, the particles, situated initially at the points A _1 fc with integer fc, move with constant velocity v. 
Then from equations (H]),© we have 

£ = NS, N = [£A] + [Xvt] (4) 

Hence 

N - [XNS] = [Xvt] 

[SXvt] SXvt 



N = N(t) = 



, jflNSX] 1-SX 
for large t. Then 



t 1-SX 



3.1 The case of nonzero drift 

Fix arbitrary e > 0, so that (X + e)S < 1. Let c e = (1- (X + e)Sy 1 (X + e)(v + e). The event A(t) = {N(t,u) > c e t} 
can be written as the union 



A(t) = U n B(n,t) 

where B(n,t) is the event that to the time t the boundary has absorbed exactly n particles, situated initially 
at the interval (0, c e t5 + (v + e)t) and more than m = c e t — n particles, situated initially at the interval 
(c e tS + (v + e)t,oo). Denote B e (t) the union of events B(n,t) with n > (A + e)(Sc e t + (v + e)t) = c e t. The 
event B e (t) has exponentially small probability, as it belongs to the event C e (t) that at time in the interval 
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(0, c e td + (v + e)t) there were not less than (A + e){5c <L t + (v + e)t) = c e t particles. The probability of the latter 
event can be estimated via Poisson distribution: 

P(C £ (t)) < C e- h W (5) 

where constant Co depends on A and e, h(e) > 0. 

Denote D c (t) the union of the events B(n, t) with n < (A + e)(5c c t + (v + e)t) = c e t. The event D e (t) belongs 
to the event that at least one of the particles, situated initially in (c e tS + (v + e)t, oo), before t will be absorbed 
by the boundary. The latter event in its turn belongs to the event J e (t), that at least one of the particles, 
situated initially in (c e tS + (v + e)£,oo), had been to the left of the point c e t8 at least once before t. Thus, 

A{t) C C € (t)U J e (t) 

The probability of the event J t {t) can be estimated from above by the expectation of the number of particles, 
situated initially in (c e i<5 + (v + e)t, oo), and visiting at least once the left side of the point c e tS before t 



P(Je(t)) < A / p(t,x)dx 

J c e tS+(v+e)t 

where p(t, x) is the probability that the particle started at x, x > c t t5 + (v + e)t, during t will visit the left side 
of the point c e t5. More exactly, 

p(t, x) = P(min(a; — vs + w(s)) < c e tS) (6) 

s<t 

From the definition of Wiener process it follows that 

P(min(a; — vs + w(s)) < c e tS) = P(mhx(—vs + w(s)) < —x + c e tS) 

s<t s<t 

= P(max(vs + w(s)) > x — c e tS) 

s<t 

To find p{t, x) we will use the well known formula for the distribution of the maximum of Wiener process with 
drift (see for example, [T], p. 935) 

P(max(« S + w(s)) > x) = 1 - $ O^^j + ( ~^^ ) ( 7 ) 

for $(z) = -±= f z e-v~l 2 dy. Then 

P(Je(t)) < A / $ dx 

J c t t6+(v+e)t \ \t J 

+ a r *»* (-*+c.t5-vt\ dx (8) 

= h(e,t) + I 2 (e,t) 

where $ = 1 — Change of variable z = x - c '*^r vt i n the first integral in the formula §8§ gives 



h(e,t) = \Vi / $(z)dz < A 

JeVi 

with ey/i > 1. 



2. 6 



Change of variables z — x c ^+ vt in the second integral in (JSJ) gives 

I 2 (e,t) = XVte- 2v2t / e 2vViz <S>(~z)dz 

J(2v+e)Vi 

= XVte' 2v2t / e 2vV ~ tz $(z)dz 

J(2v+e)Vi 
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For es/t > 1 we have 



/ j. rOO 

J 2 ( C) t) < Xj — e -(z*/2-2vzVt+2vH) 



dz 



(2w+e)Vt 



V 2tt 



Thus one can see that for e\ft > 1 



it 2 t 

= xVt<f>{eVt) < AW— e"T- 

V 2-7T 

-e"- (9) 

As the event A(t) belongs to the union of the events C e {t) and J e (t), then using the estimates §5§ and 
we get that for any sufficiently small e > 

K^^fSrHM (10) 

for some cii(e) > 0. 

Consider now the event A'(t) = {N(t,u>) < c' e t}, where c' e = (1 — (A — e)(5) _1 (A — e)(v — e). Define the event 
C' e (t), that at time inside the interval (0, c'j.5 + (v — e)t) there were not more than (A — e)(5c' e t + (v — e)t) = c' e t 
particles and let C' e (t) be the complement to C' t (t). We can write the event A'(t) as the union 

A'(t) = B' e (t)UD' e (t) 

where B' e (t) = A'(t) n C' e (t) and D' e (t) = A'(t) n C' e (t). The event B' e (t) has exponentially small probability as 
it belongs to the event C' e (t), the probability of which can be estimated via Poisson distribution: 

P(C e (t)) < C'oe-^W (11) 

where the constant Cq depends on A and e, h'(e) > 0. 

The event D' € (t) belongs to J' e (t), that at least one particle of those initially in the interval (0, c' e tS + (v — e)t) 
was not absorbed by the boundary during time t. The probability of event J' e (t) can be estimated from above 
by the expectation of the number of particles, that were initially inside the interval (0, c'^td + (v — e)t), and 
deviated more than on et its mean assignment point (if at time zero the particle is at some point x, then its 
assignment point is x — vt) at time t: 

pc e t8+(v—e)t 

P{J'M) < A / p'(t,x)dx 
Jo 

where p'(t, x) is the probability that a particle started from x, x < c? e t5 + (v — e)t, at time t will be outside the 
interval (x — vt — et, x — vt + et), that is 

i px — (v—e)t , ,o 

p (t, x) = 1 / e 2t dy 

\/2nt Jx-(v+e)t 

After substituting z = v ~^jj Vt , we find 

1 r^ 1 z 2 - 

p'{t,x) = 1 = / e~~dy = 2<$>(eVt) 

V2vr y_eVt 

Thus 

P(J'M) < C e te~^ 
Taking into account ([TT|) . we get that for any sufficiently small e > 

P ( N(t) < < X = O (e-^A (12) 



1 - {X-e)5 

for some 02(e) > 0, together with (JTUJ) this proves the first assertion of the theorem 
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3.2 Zero drift 



Here we prove the second assertion of the Theorem. 

We start with "trivial" estimate from below, we mean by this that we estimate just the number of particle 
collisions with the fixed boundary. Consider an auxiliary model where the boundary does not change its position 
while hit by particles but the particles disappear at the first collision with the boundary. Let L(t) be the number 
of collisions with the boundary for this model during time t. Evidently, we have L(t) < N(t) a. s. Denote 
Ko(t) = Ko(t,d) the number of particles, which did not reach during time t and at time t were to the left of 
the point dy/t, where d > 0. Put K(t) = N(t) - L(t). If the condition 8L(t) > d^ft holds then K(t) > K (t), 
as in the basic model all particles, situated at time t inside [0, SL(t)}, would have been absorbed by the moving 
boundary. Thus one can state that 



MK(t) > MK (t)I(L(t) > dVt) 



(13) 



where I(L(t) > d^/t) is the indicator of the corresponding event. 

Let us find the joint distribution of the variables L(t) and Ko(t). We begin with the means. Let p(t,x) be 
the probability that a particle situated initially at x, during time t will reach the boundary: 



p(t, x) = P(min(x + w (s)) < 0) 

s<t 



As 



P(min(a; + w(s)) < 0) 

s<t 



P(mmw(s) < —x) 

s<t 

P(maxw(s) > x) 

s<t 



then by formula ([7]) with v = we have 



p(t,x) = l-*(^) 



'ViJ 



2$ 



Vt) 



Let us show that for zero drift (v = 0) 



In fact, we have 



ML(t) = AW — 



ML(t) = X p(t,x)dx = 2X 



/-CO 

2XVi / $ (x) dx = X 
Jo 



Vi) 



dx 



where we integrated the latter integral by parts 



$ (x) dx 



xd§(x) 



2 ' 2 dx = 



'2tt 



Let us find MKo(t). Denote a(t,x) the probability that a particle, situated initially at x, during time t will 
not reach and moreover at time t it will be to the left of dy/t. Otherwise speaking 

a(t, x) = P(x + min w s > 0, x + Wt < dVt) — Pfmax w s < x,Wt > x — d\fi) 
[o,t] [o,t] 



The probability <j(t, x) can be written as well in the following way 

a(t, x) = Piwt > x — dy/i) — 2P(wt > x) + P(maxw s > x,Wt < x — dy/t) 

[o,t] 



G 



where due to [2], (see p. 45), we have 

/ 7y pX—dy/t pOO 

P(m&xw s > x,w t < x-dVt) = \ —~ da (2b - a)e- {2b - a)2 /2t db 

[o,t] V nt 3 J_ 00 J x 
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= $ 

Using this formula wc find 



2?rf 
X + d\/i\ 



-dy/i 

e -{2x-u)V*t du 



/+ i 5. ( x ~ d\/t\ _= / x \ , = fx + d\fb 
and then 

MK (t)= / a(t,x)dx = Vt ( / $(x)dx-2 / $<(x)dx + / $(x)<fx) = aVt 

JO \J-d JO Jd J 

where the constant 

poo poo poo pO pd 

$(x)dx - 2 / 3>(x)dx + / $>(x)dx = / $>(x)dx - / $(x)dx 



Lemma 4 Random variables L(t) and Ko(t) are independent. L(t) has Poisson distribution with parameter 
AW— , and -Kq(£) /ias Poisson distribution with parameter o~\ft. 



Proof. Let 7r„(a) be Poisson-distribution with parameter Xa, L(a,t) be the number of particles, which were 
inside (0,a) at time t = and during i were absorbed by the boundary, Ko(a,t) be the number of particles 
which were inside (0, a) at time t = 0, during time t did not reach and at time t were to the left of the point 
d\/i. Put 

q(t, x) = 1 — p(t, a;) — er(i, a;) 
Find the joint distribution of the random variables L(a,t) and Ko(a,t) : 

oo „ a „ a k rn n—k—tn 

P(L(a,t) = k,K (a,t) = m) = ^ 7r n (a)a~"^/ .../ J^P^i x it )dx H TJ <r(t, x n )dx j[ J~| q(t,x r[ )dx r[ 

n = k+m. •'° •'° i = l i = l Z=1 



where the inner sum is over all ordered arrays i\ < . . . < i& and j'i < ... < j mj of length k and m correspondingly, 
such that {«!<...< it} n {ji < . . . < j m } = 0- Here 

{n < . . . < r n - k - m } = {1, n} \ {ii < . . . < i k } U {ji < . . . < j m } 

As the inner sum consists of 



k\m\(n — k — m)\ 

equal terms (p(t,a)) k (a(t,a)) m (q(t,a)) n k m , where p(t, a) — J^p(t,x)dx, a(t,a) = Jq a(t,x)dx , q(t,a) = 
f£ q(t, x)dx, then 

oo i 

P(L(a,t) = k,K (a,t) = m) = ]T tt» (a)a~ n (p(t,a)) fc (a(t,a)) m (q(t, a))^" 1 

n—k-\-m 

{\p(t,a)) k (Xa(t,a)) m c _ Xa (Xq(t, a)) n ~ k ~ rn 

k\ ml ^—^ (n — k — m)\ 

n—k+m 

(\p(t, a)) k e' A ^' a ) (Xajt, a)) m e ~ xa ^ 
kl m\ 
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To get joint distribution of L(t) and Ko(t), let us tend a to infinity. In the limit we will get direct product of 
Poisson distributions with parameters ^\f^ and <Jy/i, as p(t,a) — > and &(t,a) —> <j\ft for a — > oo. The 

lemma is proved. 

Let us prove now that the constant in the lower bound is in fact greater than the trivial one. Choose d > 
so that d < 5ML(t)/y/t. By (fT3|) and independence of random variables L(t) and Ko(t), we have 

MK{t) > MK (t)P(L(t) > dVt) 

and 

MK (t)P(L(t) > dVt) = aVt - MK (t)P(L(t) < dVt) 
As, for a given choice of d, the probability of the event L{t) < dy/i is exponentially small 

P(L(t) < dy/i) < C 4 e-^ 

then one can assert that MKo(t)P(L(t) > d\/t) > aoVt for some constant ctq > 0, where erg < a. 

Prove now the bound from above in the second statement of the Theorem. Our goal is to estimate the 
probability P(N(t) > x\/t) from above. Take arbitrary e > 0, so that (A + e)S < 1. Subdivide the event 
A(t,N) = {N(t,uj) = N} into the union of events 

A(t,N) =U n < N B (N,n,t) 

where Bq(N, n, t) is the event that the boundary absorbed exactly n particles of the particles, situated initially 
inside (0,N6), and exactly m = N — n particles of the particles, situated initially inside (NS, oo). Denote 
Bo tt (N,t) the union of events Bo(N,n,t) with n > (A + e)SN. The event Bo^ e (N,t) has exponentially small 
probability, as this event belongs to the event Co, e (iV), that at time inside (0, N6) there was no more than 
(A + e)8N particles. The probability of the latter event can be estimated by the Poisson distribution: 

P(C .AN)) < C e~ h ^ N (14) 

where h(e) > and Co is some constant depending only on A, e. 

Denote Do }fL (N,t) the union of events Bo(N,n,t) with n < (A + e)5N. From the considerations above it 
follows that the event D e (N,t) is equivalent to the event 

m = N-n> N(l-(\ + e)6) 

At the same time the event m > N(l — (A + e)6) belongs to the event that among the particles which were at 
time to the right of the point NS, more than N(l — (A + e)S) particles during time t were at least once to the 
left of the point NS. Denote the latter event by Jo )£ (JV, t), and denote the number of particles which were 
at time to the right of the point NS but having visited the left side of NS at least once during time t (that 
is P(Jo, e (N,t)) = P(n(t) > N(l — (A + e)S)). The distribution of this random variable does not depend on 
and coincides with the distribution of the random variable L(t), defined above. By lemma [4] and the properties 
of Poisson distribution, we have the following estimate 

P(n(t) > yVt) < C ie - 0Vi (15) 

valid for some j3 > 0, y > yo > Ay^| and constant C\, depending on yo. 

As the event A(t,N) belongs to the union of the events Co, e (A) and Jo, e (N,t), then 

P(N(t)>xVi)< ]T P(C , £ (A))+ £ P(J ^N,t)) 

N>xVt N>x^/i 

By (fl4|) , for the first sum the following estimate holds 

P(C , e (JV)) = O (e"^) 

N>xVi 



8 



The second sum, by (|15|) . has the following estimate 

P(JoAN,t))= J2 P^(t)>(l-(X + e)S)N)<C 2 e'^ t 

N>xs/i N>xVt 

for some constant /3(e) > 0, x > xq — yo(l — (A + e)i5) _1 and some constant C2, depending on xq. It follows, 
that for sufficiently large x 

P(N(t) > xVt) < Cae" 7 ^ 

Hence, for t sufficiently large 

MN(t) < CVt 

The Theorem is proved. 

Proof of the second part of Lemma 1 Consider first the case of zero drift. Let us find the distribution of 
particle configuration at the time of first hitting the boundary. Put 

a(t,x,y)dy = P(x + w t 6 dy,x + min wt > 0),x > 0,y > 

s<t 

Thus, o~(t,x,y)dy is the probability that the particle which were at time t — at x > 0, during time t will not 
reach and will be in the interval (y, y + dy) at time t. As the distribution of w t coincides with the distribution 
of — w t , we have 

a(t, x, y)dy = P(x — Wt € dy, x — maxw ( > 0) = P(wt € d{x — y), maxro t < x) 

s<t s<t 

To calculate a(t,x,y) we will use the following formula for the joint distribution (see, for example [2]) 
P(w t £ da, maxti) s e db) = \ /-^r(26 - a)e" (26 " a)2/2t dad6, < 6, 6 > a 

s<t V 7rf J 

The change of variables u = (26 — x + y) 2 /2t gives 

a(t,x,y) = J^f (2b~x + y )e-^-* + ^ 2 l 2t db 

(x+ V f/2t 

e~ u du 

(x-y) 2 /2t 



L ^ e -(^-a) 2 /2t _ e -(z+y) 2 /2^ 



2vrf 

Let r be the moment of first hitting the boundary by one of the particles. 
Lemma 5 The probability density of the random variable ris 

P(t € dt) = -^e- x V^dt 

y/2^t 

Under the condition t — t the particle configuration at time t is an inhomogeneous Poisson point field with the 
rate 

poo 

\ip(t, y) = X a{t, x, y)dx = \(Q(y/Vi) - 
Jo 

where <£>(y) = ^= f*^ e~ z ' ''^ 2 dz. 

Proof. Let B±,...,Bi be an array consisting of I pairwise non intersecting intervals on the positive axis. Let 
us find the joint distribution P(r G dt,r](t, B{) — ki,i — 1,.../), where n{t, B{) is the number of particles in 
the interval B^ at time t. Let ir n (a) be the Poisson distribution with parameter Aa, r](a,t,Bi) be the number 
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of particles, which were initially inside the interval (0,a) and at time t were in Bi. Let r(a) be the first time 
moment, when one of the particles, which were in (0,a) at time t = 0, hits the boundary. Denote g(t,x) the 
probability density of the time f3 x of first hitting the point by the particle under the condition that at time 
t = this particle were at x. It is well known that (see for example [2]), that 

g(t,x) = -f=e~*' ''/ 2t 
Let us find the joint distribution of random variables 77(0, t, Bi), i = 1, I and r(a) : 

00 pa pa 

P(r(a) E dt,ri(a,t,Bi) = k u i = 1, .../) = ^ 7r„(a)a~' 1 ^/ .../ g(t,x m )dtdx m 

n=l+An+...+fei ^° J ° 

WW_dx i]s \ a(t,x h 3 ,y)dy 

3=1 s =i J v^ B i 

n — 1 — fei — ... — hi r 

TT dxr s / a(t,x Ts ,y)dy 

s=l JyeR + \u' J = 1 B 3 

where the inner sum is over all pairwise non intersecting ordered arrays m, i% t x < ... < ii,ki, *2,i < ••• < *2.fc 2 v-j 
ill < ... < ii.k t of the lengths 1, ki, k%, ...hi correspondingly. Here 

{n < ... <r n _i_ fel _..._ fc( } = {l,...,n}\({m}U{i 1 ,i < ... < i likl } U ... U < ... <ii. kl }) 

As the inner sum consists of 

n! 



ki\...kK n - l - k - ... - h)\ 

ski /»/ . r-, \ , ,1 n \\n— 1 — Aii- 



equal terms g(a, t)dt Ylj=i (&( a i ^, Bj)) 1 (a (a, t, R+ \ u' =1 £?j))™ 1 "' ' , where 

g(t,a) = / g(t,x)dx 



a(a,t,Bj) = dx a(t,x,y)dy 
Jo JyeBj 

then 

P(t(o) G dt,T](a,t,Bi) = fa,i = 1, ...I) = n n (a)a~ n 



, , fci!...A,!(n - 1 - fe x - ... - fej)' 

n— l+Ki + .-.+fc; 

I 

ff(a J t)dtn(a(a J t,B i ))^ (a(a,t 7 R + \U l J=1 B j )) n - 1 - kl -- k ' 
3=1 

,,n (Ao-(a,i,B,-)) fc ' a« 
Ag(a, J__[ , . 3!! e~ Aa x 



fe,-! 



X 



{Xa(a,t,R + \U l J=1 B 3 )y 



(n- 1- fei - - ifej)! 



j=i 3 



Taking into account that 
l 



L pa pa 

Y j a(a,t,B j ) + a(a,t,R + \U l j=1 B j )= P(fi x > t)dx = a - P(f3 x < 
• =1 Jo Jo 



t)dx 
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we have 

j pa 

-a + a(a, t, R+ \ U l j=1 Bj) = - V &(a, t, Bj) - / P(/3 X < t)dx 

Thus 

P(r(a) e dt, V (a,t,B t ) = k i: i = = Xg(a, t)e~ SS ne*<t)d* dt J| (Afr(a, t, B 3 -))*' e -A» ( ,, t , B ,) 

To obtain the joint distribution of r and ry(t, £?,), we tend a to infinity. Then 

#(a,f) — > ^ g(t,x)dx = \J^ t 

j P((3 x <t)dx ^ P{p x <t)dx = sj^ 

a(a,t,Bj) — > a(t,Bj) = / c?y / a(t,x,y)dx 

JyeBj JO 

= / cW(*> y) 

where tp(t,y) = $(y/Vi) - ^(-y/y/t). Hence 

P(r e dt, »?(*, SO = fci, i = 1, ...I) = -^=e" V?df (^M!i e -Aa(t,B J ) 

\/2,7rt .. 

The lemma is proved. 

Let us prove that the random variable fci can take infinite values with positive probability. Assume that 
t = t. Denote mi the number of particles in the interval [0, 8} at time r = t, mi the number of particles at time 
r = t in the interval D\ = [8, (mi + 1)5], if mi > 0, and denote the number of particles at time r = t in the 
interval D 2 — [(mi + 1)8, (mi + 1)8 + mi8], if mi > 0, m 2 > 0, etc. The sequence of pairs [n = (mj_i, mj) is a 
discrete time homogeneous Markov chain. Transition probabilities are 



(A*(t,£>or'+i r _ wu?< , 

m i+ i! 



i 



The states (m, 0) are absorbing. Random variable ki is defined by the random stopping time of the Markov 
chain fi^. The boundary stops, when the chain hits some of the absorbing states. We introduce Lyapounov 
function f(fii) = m^. Its one step increment equals 

M(f(fj, i+ i) - f(m)\ni = (mi-i,mi)) = Xa(t,D l )-m l 

/• (mi _i+l) 5 +rrti 5 

= A / ip(t,y)dy - mi 

J{mi--L+V)5 

As tp(t, y) — > 1 for y — » oo and XS > 1, then there exists such e > 0, that for sufficiently large mj_i 

M (f(Hl+l) - = K-b™.)) > £ ™J 

This is the transience of the Markov chain pn. Thus, P(ki — oo) > 0. 

The case of nonzero drift is considered similarly. Moreover, as we consider the lower bound, it is intuitively 
clear that the drift can only increase the probability that the boundary explodes to infinity. 
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